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We derive the spherical field formalism for fermions. We find that the spher- 
ical field method is free from certain difficulties which complicate lattice 
calculations, such as fermion doubling, missing axial anomalies, and compu- 
tational problems regarding internal fermion loops. 
[PACS numbers: ll.10.Kk, 11.15Tk] 



1 Overview 

Spherical field theory is a new non-perturbative method for studying quan- 
tum field theory. It was introduced in [1] and was used to describe the 
interactions of scalar boson fields. In this paper we show how to extend the 
spherical field method to fermionic systems. 

The central idea of spherical field theory is to treat a (^-dimensional system 
as a set of coupled one-dimensional systems. This is done by expanding field 
configurations of the functional integral in terms of spherical partial waves. 
Regarding each partial wave as a distinct field in a new one-dimensional the- 
ory, we interpret the functional integral as a time-evolution equation, where 
the radial distance in the original theory serves as the time parameter. For a 
purely bosonic system the time-evolution equation corresponds with a mul- 
tidimensional partial differential equation. In the case of a purely fermionic 
system, we find that the time evolution is described by a system of first-order 
ordinary differential equations. In future work we will study mixed systems 
with both bosons and fermions which are described by coupled partial dif- 
ferential equations. 

Unlike lattice methods, spherical field theory yields an expansion which, 
at any order, corresponds with a continuous system. It is therefore able to 
avoid problems associated with discrete approximation methods.^ There is 
no doubling of fermion states, and we find the correct axial anomaly. Fur- 
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thermore internal fermion loops present no special computational difficulties 
and is included in the dynamics of the time-evolution equation. Detailed 
examples of such calculations will be presented in a forthcoming paper. We 
anticipate that spherical field methods will be useful in the study of non- 
perturbative fermionic systems, especially chiral fermions and phenomena 
related to fermion loop processes. 

The organization of this paper is as follows. We begin with a brief descrip- 
tion of Grassmanian path integrals and the fermionic analog of the Feynman- 
Kac formula. We then generate the spherical expansion for free fermion the- 
ory in two Euclidean dimensions with sources and derive the spherical field 
time-evolution operator and generating functional. By functional differen- 
tiation with respect to the sources, we obtain the spherical field formalism 
for general interacting theories. Next we check that spherical field theory 
produces the correct axial anomaly. We then show how to write the time- 
evolution equation as a matrix system and comment on the utility of spherical 
field methods in studying fermionic systems. Although our analysis is done 
in two dimensions, the extension to higher dimensions is straightforward. 



2 Grassmanian path integrals 

Let V>j(i) and ipj(t) be Grassman-valued functions where i,j — l,---,N. Let 
Vfyi, ipj,t) be a polynomial in ^ and ipj, ordered such that all ^'s are placed 
on the right and all ?/>'s are placed on the left. In [[J it is shown that 
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This is the fermionic version of the Feynman-Kac formula. A more recent 
derivation using fermionic coherent states can be found in || and [|J . The an- 
tiperiodic boundary conditions imposed at tj and tp follow as a consequence 
of computing the trace. We note that these are in fact special conditions. 
More general boundary constraints produce ambiguities which depend on the 
specific discrete approximation used to obtain the continuum limit. 

It is not clear that such antiperiodic boundary conditions can be general- 
ized in a coordinate-independent manner for functional integrals over higher 
dimensional regions. The rigorous theory of Grassmanian functional inte- 
gration has not developed to the point where we can answer such questions. 
Nevertheless functional integration is a convenient method for deriving useful 
field-theoretic results, although in a somewhat heuristic fashion. In this anal- 
ysis we use the functional integral to deduce the spherical field formalism for 
fermions. Although we will be careless with regard to boundary conditions, 
in the end we explicitly check that the spherical field method produces the 
correct generating functional for free field theory. By functional differentia- 
tion with respect to the external sources, we conclude that the spherical field 
formalism is valid for general interacting theories. 

3 Spherical fermions 

Let us consider Euclidean field theory in two dimensions. We will use both 
cartesian and polar coordinates, 

t— (t cos 9, t sin 9) = (x,y). (4) 

In Euclidean space the gamma matrices satisfy 

{y,y} = -2*« (5) 

and we choose the representation 

7 = io. (6) 

Let us start by constructing the spherical field Hamiltonian. We first decom- 
pose the fermion fields, 
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n=0,±l,- 
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a- V = 

we have 
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The Euclidean action for free field theory with external sources, r\ and fj, is 
S = -i J dQdtt ($(i7 • V - m)^ + ^ + #) . (11) 
In terms of partial waves, 4 



= -% f dtt 
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The generating functional is therefore 

J j exp 1 jf" (ft t G n J , 

where G n is defined as 



(12) 

(13) 
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If we now define ^ = t^ l n , the generating functional is, up to an overall 
constant, 

/ (n^Wn) ^p | dt _E , (15) 
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where G' n is 

+ - ^-i (# - (i6) 

Our goal is to find an equivalent expression for (|I3|), in analogy with the 
Feynman-Kac formula ([1]). We start by defining a linear vector space. For 
each finite subset 

' n = 0, ±1, ±2, • - ■ 1 



5 C 



(17) 



we assign a vector |5) satisfying the following orthogonality and normaliza- 
tion conditions, 

(m = [ I ssts>. (18) 

For later convenience we define a lexicographic order, namely, 

4 < 4 (19) 

if and only if 

n < n (20) 
or n = n , i = I , and i' = \ . 

Let £ be the linear space spanned by all such vectors \S) . Let us define 
operators and a % ~ by the following relations, 

<;][S) = { (-l)#|°S-<> if s£ el (21) 



< + \S) 



if 4 e 5 
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where # is the number of elements in 5* which are less than s l n . Comparing 
(15) with ([!]), we make the correspondences^ 

V'l^i'n-i - (22) 



5 We denote the conjugate of a l ~ as a^ + . Although corresponds with a partial wave 
with orbital angular momentum n, we note that corresponds with a partial wave with 
orbital angular momentum — n — 1 or — n + 1, depending on i. 
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and define 
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The traces in ([£3) are performed over the space E. In practise it is not 
necessary to explicitly compute these traces, since only the ground state 
projection contributes. For m — 0, however, the ground state is degenerate 
and the e — > prescription picks out the correct ground state. We will refer 
to (^) as the fermionic spherical field ansatz. For notational ease we will 
suppress the e terms. We now show that Z[fj,r]] is the correct generating 
functional for free field theory. 

We note that S can be decomposed as a tensor product space with the 
identification 



\S) 
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Since H n (r),r),t) acts upon only the nth component in the tensor product, 
we can write Z[fj,rj\ 9S 9 product 
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In (27) the trace is performed over the four- dimensional space spanned by 
the vectors 

10), |Rh}), |{4>), |{44 + i}}- (28) 

In Appendix 1 we derive the result, 
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where the matrix M n (t 1; t 2 ) is defined asf] 

mK n (\m\ ti)I n {\m\ t 2 ) \m\ K n (\m\ ti)I n+1 (\m\ t 2 ) 
m\ K n+1 (\m\ ti)I n (\m\ t 2 ) mK n+x (\m\ t x )I n+x (\m\ t 2 ) 

ml n (\m\ t x )K n (\m\ t 2 ) - \m\ I n (\m\ t x )K n+1 (\m\ t 2 
- \m\ I n+1 (\m\ t x )K n {\m\ t 2 ) ml n+1 (\m\ t x )K n+x (\m\ t 2 ) 
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for m — 0. 

Let us now compare these results with the known results for free field 
theory. The two-point free field correlator is 
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From (33) and (|3~4] ) it is straightforward to recover the spherical correlation 
functions in (|30| ) and (|31~D . We conclude that the fermionic spherical field 
ansatz produces the correct generating functional for free fermions. These 
correlation functions are part of the spherical Feynman rules for fermions. 
For future reference we have written these in a more convenient format in 
Appendix 2. By functional differentiation with respect to the sources rj and fj, 
we conclude that the spherical field formalism is valid for general interacting 
theories. 



e Ii and Ki are the i th order modified Bessel functions of the first and second kinds 
respectively. 
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4 Axial anomaly 



We now show that the spherical field formalism yields the correct form for the 
axial anomaly. We consider free massless fermions, again in two Euclidean 
dimensions. Let us define 



S^(t) = (0 V^i)A u (0) 







(35) 



where and are the vector and axial vector currents, and the subscript 
E is intended as a reminder that our Euclidean correlation function is defined 
as the analytic continuation of the corresponding time-ordered function in 
Minkowski space. 7 We note that 
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The one-loop process corresponding with S^ u carries a logarithmic divergence 
proportional to e^ v , and the regulated value of S^ v will depend on our def- 
initions of V 11 and A^ as operator products. In our discussion here we will 
remove this ambiguity by considering the symmetric combination S^ v + S UfJ- . 
Let us define vector and axial vector currents, 



Expanding the currents in terms of partial waves, we have 
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From (39) and (40), we have 
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Integrating with any smooth test function we recognize these terms as deriva- 
tives of the Dirac delta distribution, 
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We conclude that 
J d 2 te l ^[(o\v^(?)d u A u {0) 
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\e^ Vv . (47) 



If we now choose to maintain conservation of the vector current, we have 



J d 2 te l ^(o\v»{i)d u A u {0) 







(48) 



which is the desired result for the axial anomaly (see [6]). This should not 
be surprising. With d v A v placed at the origin, the calculation we have done 
is the same as that of standard field theory in position space. For the case 
when d v A v is not at the origin, our spherical field calculation generates an 
expansion of 



S^(t- i>) = (0 \V»{t)A v (?) 







E 
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in terms of sums of spherical waves. The partial sums of this expansion 
converge pointwise (except at t = i 1 ) to the result (44), and we again get the 
correct axial anomaly. 
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5 Matrix representation 



Spherical fermion fields can be studied using ordinary matrices. Unlike the 
spherical bosonic system which corresponds with a multidimensional partial 
differential equation, the spherical fermionic system corresponds with a set of 
coupled first-order ordinary differential equations. We illustrate some basic 
methods here using the free fermion system. 
Let us define the following column vectors 
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We now write the Hamiltonian and creation and annihilation operators as 
matrices, 
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As an illustrative example we calculate the correlation function 

<0|W(*)$(0)|0) B (53) 

for the case m^fl. Let 

U(t 2 ,fi) = Texp j-^ 2 dt# (0,0,t)j. (54) 
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Using the fact that 
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which agrees with (|3T|). 

Calculations in interacting theories are done in a similar manner. There 
we encounter new terms in the Hamiltonian which can be found by functional 
differentiation with respect to the sources rj and fj. Several detailed examples 
of interacting systems will be presented in later work. 

6 Summary 

In this paper we have extended the spherical field formalism to include 
fermionic systems. Since spherical field theory deals with continuous sys- 
tems, it avoids problems associated with discrete approximation methods, 
such as fermion doubling and cancelled axial anomalies. The spherical field 
method should therefore be useful in studying chiral fermion systems. 

We have shown that the time evolution of spherical fermion fields can 
be modelled using matrices and is described by a set of coupled first-order 
ordinary differential equations. We recall that in lattice field theory Grass- 
manian variables are associated with each fermionic degree of freedom at 
each lattice site. In spherical field theory, however, Grassmanian creation 
and annihilation operators are associated with each spherical partial wave. 
Although we have not analyzed interacting systems in this work, we antici- 
pate that for comparable levels of accuracy spherical field theory will require 
manipulating much smaller anticommuting algebras. Since large sets of anti- 
commuting variables present serious computational problems, spherical field 
theory may yield significant improvements in the numerical calculation of 
fermionic interactions. 

Acknowledgements 

The author is grateful to Howard Georgi for useful discussions regarding 
this work and Robert Shrock for providing the impetus and encouragement 
to investigate fermions using spherical field theory. The author also thanks 
the Harvard physics department, where this work was started. 

Appendix 1 

In this appendix we derive the results (29), (|30D, and fl3T|). Let us define 

W n = \nZ n . (63) 
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Using the notation, 



we note that 
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Differentiating with respect to t we find 
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Differentiating again we get 
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(71) 



13 



For m^O the general solution for W n is 

J dhdt 2 [ hfjlnih) hfj^ih) ] M n (t 1? t 2 

J dttK n (\m\t)fjl n (t),---,J dttl n+1 {\m\t) v l n+1 (t) 
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where F is the general homogeneous solution and M. n (ti,t 2 ) is given by 
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for m — 0. The function F depends on eight anticommuting variables and is 
therefore a polynomial of degree at most eight. It is straightforward to check 
that for ti 7^ 0, the limits 



lim - f - , Z n -A— 



r)=7]=0 



and lim 



a 17 



r)=r)=0 



(75) 



(76) 



are well-defined. Corresponding limits for the other two-point correlators are 
also well-defined, and the same holds true for other correlators such as 
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(78) 



provided that t 1 ,t 2 ,t 3 7^ 0. From this and the fact that W n vanishes when 
rj = fj = 0, we conclude that F = and we obtain (29). 
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Appendix 2 

The following is a list of the spherical correlation functions. For m 7^ 



we 
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